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Abstract 



^ ! One random spin-| XY chain that after Jordan- Wigner fermionization 

Q \ reduces to the extended Lloyd's model is considered. The random-averaged 

^ \ one-fermion Green functions have been calculated exactly that yields thermo- 

f~^ ' dynamics of the spin model. 
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An idea to exploit Lloyd's model |^ for examining the tliermodynamical properties 
of random spin-^ XY chains belongs to H.Nishimori [Q. He noted that after Jordan- 
Wigner trick the Hamiltonian of isotropic XY model with random lorentzian transverse 
field describes tight-binding spinless fermions with diagonal lorentzian disorder. Since 
the random-averaged one-fermion Green functions for such model were found exactly by 
P.Lloyd, one can obtain the thermodynamics of random spin system via the averaged 
density of states. Later the treatment presented in [Q was generalized for the cases of 
alternating bonds [Q and additional intersite Dzyaloshinskii-Moriya interaction 0. 

On the other hand, W.John and J.Schreiber suggested an extension of Lloyd's method 
to off-diagonal disorder that was successfully used in the study of disordered systems 
[7-11]. The idea of the present communication is to exploit Lloyd's model with off-diagonal 
disorder for analysis of thermodynamics of the corresponding random spin-^ XY chain. 
Similarly to we were able to calculate exactly various thermodynamical quantities, 
although found somewhat different results of influence of randomness on these functions. 

We consider N spins | arranged in a circle with the Hamiltonian 

TV N 

n=l n=l 

where Qn is a transverse field at site n and J„ is the interaction between the sites n and 
n -|- 1. The latter are taken to be random with a probability distribution density 

TV ^ Y 

that is the product of lorentzian distribution densities at sites that are centered at Jq with 
the width F. In order to treat the model (1), (2) in exact manner the transverse field f2„ 



at each site must depend on surrounding intersite interactions in the following way 

r2„ — fio = fl ( ~ — 7. 1 — —^^ ) ; fl is real, | a |> 1, (3) 

where Qq is the averaged transverse field at site. 

Really, by Jordan- Wigner transformation from operators s • = sj ± isj to Fermi op- 
erators Cj, Cj the Hamiltonian (1) becomes 

H = R- + BP+, 

1 TV N N J 

n = —— 2_^ \ln + 2^ ^t-nCnCn + 2^ — (^C„ C„,+l — CnC^^ij , Cn+N = Cn, C^+jy = C„ , 
^ n=l n=l n=l ^ 

1 _|_ p N 

B^-Jr,{4c,-c^ct), P^^-—, P=n(-2<)-(4) 

^ n=l 

For calculation of thermodynamical properties of the model (1) one can omit the boundary 
term B [Q, and hence one faces with one- dimensional version of Anderson's model with 
the off-diagonal disorder considered by W.John and J.Schreiber. 

In order to study thermodynamics one should diagonalize the bilinear in Fermi op- 
erators form H' (4) by canonical transformation rjk = J2n=i QknCn with real gun that 
satisfy the equations kkgun = J2f=i9kiAin with Aj = ^Aj + l^i^i.i+i + |Ji-i5j,i-i, and 
Y.^=i9kigpi = 4p, Y.p=i9pi9pj = ^ij obtaining in result H' = Y.p=i^piVpVp - |)- The 
density of states p{E) = -^ Y.p=i ^{E — Ap) determines thermodynamics for certain re- 
alization of random intersite interactions. For example, the Helmholtz free energy per 
site is given by / = —\^dEp{E)\n.{2ch.^). The Helmholtz free energy averaged over 
random realizations is given by the same formula only with the random-averaged density 



of states p{E), where the averaging is defined by (...) = J dJi...dJNp{Ji, ..., Jn){---)- 

On the other hand, the temperature double-time Green functions TJ^^t) = ^i6{±t) < 
{Vpit),vt} >> ^W = ^I-oodEe-'^rj^iE ± le), e -> +0 yield the density of states 



for a certain random realization: p{E) = -^ Y,p=i -f^^^'^^JpiE ± ie) . p{E) can be 
rewritten in terms of Green functions G^j^(t) = ^i9{±t) < {c„(t),c^} > as p{E) = 
T^jjYljLi'^'^GjjiE±ie) since T^g{t) = YliLiJ2f=i9pi9qjGfj{t). In result the averaged 



density of states is determined by the averaged Green functions G^^{E) via the relation 



p{E) = T^lmGUE) 



Finally, following one can derive the exact expression for G^^{E). First it is 
necessary to write a set of equations for G^j^{E±ie) that follows from equations of motion 
for G'^^it) and then to average these equations using contour integration in complex 
planes of J„s. Under the imposed condition (3) on the basis of Gershgorin criterion one 
can state that for a > 1 the retarded (advanced) Green function cannot have a pole in 
lower (upper) half-planes of J„s, whereas for a < — 1 in upper (lower) half-planes of J„. 
Therefore, every contour of integration should be closed in the half-plane where there 
is only the pole originated from lorentzian distribution, and after trivial use of residues 
one gets a set of equations for the averaged Green functions that possess already the 
translational symmetry and hence may be solved in a strandard way. The final result for 
the averaged Green functions reads 

\n~m\ 



Wx'^—y'^—x 



Gtm{E) = ^ ^=^= (5) 

with X = E — VlQ±i\ a\V, y = Jq^ isgn(a)r. 

The obtained averaged Green functions (5) permit to study thermodynamics of spin 
model (l)-(3). Really, the required averaged density of states that follows from (5) reads 

1, 1 



p{E) = =F— Im 



^ ■l{E-Vto±i\a\Tf -{J^T isgn{a)Tf 



A={E-n^f + {l-\a\^)T''-Jl B = 2T[\a\{E-n^)+sgn{a)J^]. (6) 

The entropy and specific heat can be calculated by formulae 

(7) 



dEp{E) 
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dEp{E) 2 



Vch^ 



f3E 



Due to the noteworthy property of (6) -^p{E) = —-^p{E) one can express transverse 
magnetization and static transverse linear susceptibility through the density of states 



^=<^E^^> = -^/ dEp{E) th^, (9) 

dmZ „ f .„ , „, 1 



Let us discuss the obtained results. In the absence of randomness (F = 0) (6) reduces to 



the well-known result: piE) = -—r^, — = if \ E — Qr, \<\ Jn \ and p(E) = otherwise. 

The isotropic XY model in random lorentzian transverse field treated by H.Nishimori 
may be obtained in the limit F ^ 0, | a | F = const = Fj^. The model in question (l)-(3) 
essentially differs from that model: the density of states (6) in contrast to the case of 
diagonal disorder is not symmetric with respect to the change E — VLq -^ —(E — Qq). 
However, it remains the same after the replacement E — Qq -^ —{E — fio), a -^ —a, or 
E — VIq ^> —{E — Qq), Jq -^ —Jo, since the simultaneous change of signs of Jq and a in (6) 



does not affect p{E). For convenience hereafter will be put Jq = 1. The above-mentioned 
symmetry of the density of states can be seen in Fig.l, where the averaged density of 
states (6) for F = 1 is displayed. The density of states for non-random case is depicted in 



Fig.l by dashed lines. For large | a \ due to disorder the edges of the zone are completely 
smeared out; for | a |~ 1 the disorder results in smearing out mainly of one edge of the 



zone. Some consequences induced by this dependence of p{E) on a for F 7^ will be seen 
in the behaviour of thermodynamical quantities. 

The results of numerical calculations of thermodynamical quantities for F = 1 and 
few values of a are presented in Figs. 2-5, namely, the temperature dependences of entropy 
(7) (Fig.2), specific heat (8) (Fig. 3) and static transverse linear susceptibility (10) (Fig. 5) 
and the dependence on averaged transverse field at low temperatures of the transverse 
magnetization (9) (Fig. 4); the curves that correspond to non-random case are depicted in 
these figures by dashed lines. The influence of randomness on thermodynamics is mainly 
rather typical. It leads to weak deformation of the curve entropy versus temperature with 
decreasing of entropy at high temperatures (Fig.2), broadening and decreasing of the peak 
in dependence specific heat versus temperature (Fig. 3), smearing out of the cast in the 
mj versus Qq curve at T = for Qq = Jq and nonsaturated transverse magnetization 
at any finite transverse field (Fig.4), suppressing of static transverse linear susceptibility 
versus temperature curve (Fig. 5). However, as can be seen in Figs. 2-5 the infiuence 
of disorder, especially for small a, essentially depends on the sign of a. Particularly 



interesting is the case of strong asymmetry in the density of states p{E) when | a |^ 1. 
From mathematical point of view the dependence of computed quantities on temperature 
and averaged transverse field and the well-pronounced difference between the cases a ~ — 1 
and a K. \ can be understood while bear in mind that these quantities according to (7)- 



(10) are the integrals over E of the products of p{E) depicted in Fig.l by the functions 
with evident dependence on E at different 13. It is interesting to note that for some 



Hamiltonian parameters and temperatures even the large randomness (controlled by F) 
almost does not affect the observable thermodynamical quantities. This can be nicely 
seen in Figs. 2-5. 



It is worth to underline that the asymmetry of p{E) leads to the appearance of nonzero 
transverse magnetization mj at zero averaged transverse field Qq. As it can be seen from 
(9) ml = at T = 0, fio = if /°^ dEp{E) = J^ dEp{E). This is evidently true for a 



symmetric density of states p{E) (as in the case considered by H.Nishimori) but is not 
obvious in the case in question (6). The difference between the integrals J_^ dEp{E) and 



/o°° dEp{E) can be clearly demonstrated by numerical finite-chain calculations [0 as a 
difference between the numbers of negative and positive eigenvalues oi N x N matrix 
1 1 Aij 1 1 Ap, denoted by A/1 and JV^ respectively, for certain realization of random model 
(l)-(3). For a realization of random chain (l)-(3) of 1000 spins with Qq = 0, Jq = 1, T = 1 
that gives ^ E^=i Jn = 1-009 we found that for a = -5 A/l = 495, A/V = 505, for a = -2 
A/"- = 470, A/V = 530, for a = -1.01 A/1 = 402, A/V = 598. Another random realization 
of this chain with ^ E^=i Jn = 0.986 yields for a = -5 Af- = 503, Af+ = 497, for a = -2 
U- = 471, A/V = 529, for a = -1.01 A/1 = 408, A/V = 592. The transverse magnetization 
for certain realization at T = is given by rn^ = ~~^ + and one finds a good agreement 
of calculated in such a manner — mj with the results depicted in Fig. 4. 

To summarize, this paper is devoted to thermodynamics of spin-| isotropic XY chain 
with random lorentzian intersite interaction and transverse field that depends linearly on 
the surrounding intersite interactions (l)-(3). The derived exact expressions for the aver- 
aged density of states (6) and thermodynamical quantities (7)-(10) seems to be interesting 
from academic point of view since they permit to understand the disorder effects and from 



applied point of view since they may be used as a testing ground for approximate methods 
of spin systems with off-diagonal disorder. 

Unfortunately, the obtained results do not permit to calculate exactly the averaged 
spin correlation functions because such calculation requires the knowledge of averaged 
many-particle fermion Green functions. Spin correlations and their dynamics may be 



examined using exact finite-chain calculations developed in ||1J, [15 . 
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List of figure captions 



Fig.l. The averaged density of states (6) p{E) vs. E — Qq- 

Fig. 2. Tlie entropy s (7) vs. temperature 4. 

Fig. 3. The specific heat c (8) vs. temperature 4. 

Fig. 4. The transverse magnetization —m^ (9) vs. transverse field Qq at low temperature 
(1 = 0.001). 

Fig. 5. The static transverse linear susceptibility —Xzz (10) vs. temperature 4 at r2o = 0.5. 



